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_  £  0 

This  paper  considers  solutions  of  0  =  e  \  a.  . (x)u  + 

•  •  1  «  I  A  •  A  • 

1,J  1  J 


0  0  2 
b.&Otu^in  a  bounded  domain  which 

•  1  1  ' l:  -  ■  -- 


is  bounded 


in  e  >  0.  We  assume  that^  ae  ,  be  b^ — that  all 

solutions  of  the  ODE  converge  to  a  single 

*  ,V.C.  <1  t 

linearly  asymptotically  stable  critical  point  in  0,  without 

leaving  ft),  giva  a  jaroof^  based  on  the  standard  probabilistic 

interpretation  of  uE,  of  an  exponential  leveling  property;-'4" 

sup  |u  (x)  -  u  (y) |  <  e’  '  for  some  5  >  0  which  depends  on  the 
'x,y£K 

compact  set  K  c  fi. 
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EXPONENTIAL  LEVELING  FOR  STOCHASTICALLY 
PERTURBED  DYNAMICAL  SYSTEMS 

I :  Introduction 

Consider  a  deterministic  system  described  by  an  ordinary 
differential  equation  in  ]R^: 

(1.1)  dx°(t)  =  b(x°(t))dt. 

A  natural  model  for  the  behavior  of  this  system,  when  subjected 
to  a  small  stochastic  perturbation,  is  the  diffusion  process 
described  by  the  Ito  equation 

(1.2)  dx£(t)  =  b(xe(t))dt  +  /e  o(xE(t))du>t, 

“>t  being  a  Brownian  motion  in  IR^.  Applications  of  this  type  of 

model  can  be  found  in  Ludwig  [6],  Schuss  [10]  and  Matkowsky  and 

Schuss  [9].  Several  aspects  of  the  asymptotic  behavior  of  x£(*) 

as  e  -*■  0  are  of  interest.  Consider  in  particular  a  bounded 

£ 

domain  ft.  If  denotes  the  exit  time  of  x  (•)  from  ft  and 

£  £ 

Ex  the  expectation  for  the  solution  of  (1.2)  subject  to  x  (0)  = 

then  (under  appropriate  regularity  assumptions) 

uG(x)  =  E£[f(x£(in))] 
is  the  solution  of  the  Dirichlet  problem 
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(1.3)  0  =izfe[u]  =  |  l  a..(x)u  +  l  b.(x)u  in  ft 

L  -  j  13  .  i  X;i 

with  u|  =  £. 

1  3ft 

T  E 

(Here  a  =  ao  ) .  The  behavior  of  u  as  e  +  0  depends,  of 
course,  on  the  nature  of  the  trajectories  of  (1.1)  which  start  in 
ft.  One  of  the  more  interesting  cases  is  when  all  deterministic 
trajectories  starting  in  ft  remain  in  ft  and  approach  a  unique 
stable  point,  at  the  origin  say.  Because  all  continuous  solutions 
of  the  reduced  equation 

0  =  l  b.  (x)ujj 

i  l  x£ 

in  ft  are  constant,  one  expects  that  ue  approaches  a  constant 
function,  or  at  least  somehow  "levels  out".  We  prove  here,  under 
modest  assumptions,  that  this  leveling  does  occur  and  at  an  ex¬ 
ponential  rate: 

sup  jue(x)  -  ue(y)  |  <  e~6//e 
x.yeK 

for  any  compact  K  c  ft,  some  6  >  0  and  all  sufficiently  small  e. 

In  many  cases  much  more  is  known.  Matkowsky  and  Schuss  [8] 
presented  a  formal  calculation  to  show  that  ue  converges  to  a 
constant  function  and  derived  a  formula  for  what  this  constant  should 
be.  Kamin  [5]  and  Devinatz  and  Friedman  [1]  gave  rigorous  proofs 
of  this  in  cases  where  jzfe  has  a  self-adjoint  form 


3 


(1.4) 


J/e[u]  =  e*^e 


9x . 


,  ^/e 

Ce  aijV 


In  [4]  Kamin  showed  that  the  formal  calculation  of  Matkowsky  and 

Schuss  for  (1.3)  is  correct  provided  the  solutions  of  certain 

auxiliary  first  order  PDE's  exist  and  are  sufficiently  smooth. 

The  fundamental  work  of  Ventcel  and  Freidlin  [12]  also  establishes 
0 

that  u  converges  to  a  constant  for  (1.3)  in  the  case  that 
the  variational  distance  V(0,y)  which  is  central  to  their 
treatment  attains  its  minimum  over  y  €  3ft  at  a  unique  place. 

Actually  the  in  (1.4)  is  of  a  more  general  form  than 

(1.3): 


(1.5)  i/£[u]  =  |  l  a..u  +  l  bCu 

e  0  e 

with  b  -*■  b  as  e  -*■  0.  In  this  context  the  solutions  u  may 

not  converge  to  a  constant.  Indeed,  in  [1]  the  authors  presented 

the  two  examples 


(1.6a)  e(x+2)u"  -  x(x+2)u'  =  0 

(1.6b)  e(x+2)u"  +  (e-x(x+2))u’  =  0 

on  [-1,1],  both  with  u(-l)  =  0,  u(l)  =  1.  They  observed  that 

0  1  0  T 

u  +  j  for  (1.6a)  and  u  -*■  j  for  (1.6b).  If  we  combine  these 
two  examples  as 
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(1.6c)  e(x+2)u"  +  (£  Sin(-gO  -  x(x+2))u'  =  0, 

we  get  an  example  of  the  type  (1.5)  for  which  uE  does  not  converge. 
The  result  proved  here  is  for  of  the  form 


[2]  =  £  a  .e  .u  +  £  bfu 

1  2  i*  i,3  x.x.  i  x. 


The  a8,!)6  are  required  to  converge  to  a^,b°  as  e  -t-  0 .  This 
£ 

form  of  encompasses  all  the  cases  (1.3)- (1.6)  mentioned. 

£  £ 

The  boundary  function  u  =  f  is  allowed  to  be  e-dependent 

1  an 

and  is  required  only  to  be  bounded  (in  both  x  and  e)  and 
measurable  (Borel) ,  but  need  not  converge  with  e. 

Section  2  contains  the  technical  assumptions  and  the  statement 
of  the  main  theorem.  Sections  3  and  4  are  devoted  to  a  bound  on 
hitting  probabilities  which  is  the  cornerstone  of  our  proof.  The 
proof  of  the  theorem  is  given  in  Section  4  also.  Section  5  contains 
two  additional  remarks. 
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II :  Technical  Assumptions  and  Statement  of  the  Main  Result 
The  domain  ft  c  is  assumed  to  be  bounded.  To  treat 

£ 

u  as  the  solution  to  an  elliptic  boundary  value  problem  with 
£  £ 

u  I  =  f  ,  a  specified  continuous  function,  one  might  also  want 
1  3  a 

2 

to  impose  the  requirement  that  3ft  be  C  .  The  probabilistic 

e 

definition  (2.1)  of  u  renders  this  unnecessary  however.  The 
assumptions  on  the  coefficients  are  as  follows: 
e  0 

a)  a  (x) ,  a  (x)  are  Lipschitz  (or  just  Holder)  continuous  in  x 
uniformly  in  £,  positive  definite  symmetric  d  x  d 
matrices  on  ft  and  |  ai .  -  a^  |  •*  0  uniformly  on  ft 

as  e  F  0 ; 

b)  be(x)  and  b°(x)  are  in  (^(ft),  |bE-b°|  and 

£  f) 

| b  -b  j  (i  =  l,...,d)  all  converge  to  0  uniformly 
xi  xi 

on  ft  as  £10; 

c)  for  any  solution  to  x ^  (t)  =  b^(x^(t))  with  x^(0)  £  ft, 

x^(t)  £  ft  for  all  t  >  0  and  lim  x^(t)  =  0; 

0* 

9bV(0) 

d)  the  matrix  B  =  ( — — ]  is  stable,  i.e.  all  its 

5 

eigenvalues  have  negative  real  parts. 

For  a  specified  x  €  ft,  xE(t)  is  a  Markov  diffusion  process  with 
£ 

x  (0)  =  x  and  differential  generator 


& 


e  r.  e 

=  t  l 


2  i,j  3x.3x. 


+  I  b 

i 


i  JxJ- 


£  £ 

For  definiteness,  one  can  think  of  a  ,b  as  being  extended  to 
all  of  1R^,  the  process  xe(*)  then  being  obtained  as  below  for 
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all  t  <  +“.  We  are  only  concerned  with  x£(t)  for  t  5  tq»'  ^ 
however,  which  does  not  depend  on  this  extention.  If  one  likes, 
x  (t)  can  be  considered  as  the  solution  to  a  stochastic  differ¬ 
ential  equation  (see  [7]  or  [11]) 

dx£(t)  =  b£(x£(t))dt  +  /e  ae(x£(t))dwt,  x£(0)  =  x 
6  £  £  T  £ 

if  a  =  0  (0  )  where  0  is  Lipschitz.  Alternately, 

0 

x  (t)  can  be  discussed  directly  via  the  martingale  problem 

associated  with  ;  [11].  (Continuity  of  coefficients  is 

sufficient  for  that  treatment . ) 

£ 

The  boundary  functions  f  (x)  are  assumed  to  be  bounded  in 

e 

x  and  e  >  0  and  measurable  on  9ft.  The  u  (x)  are  now  defined 
by 

(2.1)  u£(x)  =  E£[f£(x£(Tn))]. 

£  2 

It  can  be  shown  that  u  £  C  (ft)  and  satisfies 

i?f£[u£]  =0  in  ft. 

_  £ 

Indeed,  on  any  ball  B  with  B  5  ft  it  is  true  that  u  is  the 

0  1 

Perron  solution  corresponding  to  the  boundary  data  u  .  Since 

2  9B 
Perron  solutions  are  c  on  the  interior  of  their  domains,  for 

bounded  measurable  data,  ([3],  Theorem  6.11)  it  follows  that 

u  €  C  (ft).  The  boundary  behavior  of  u  does  not  concern  us; 

only  the  boundedness  in  e,x  is  necessary  for  our  arguments  below. 
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There  is  one  more  condition  that  we  will  need  when  proving 
Theorem  2  below.  In  deriving  (4.6)  we  will  use 

e  0 

(2.2)  — — (x  j  =  o(l)  as  £4-0  uniformly  in  ft. 

We  know  that  b^(0)  =  0,  so  the  above  will  follow  from  the  con- 

£  0  E 

vergence  of  b  to  b  if  the  further  condition  b  (0)  =  0 

xi  xi 

is  true.  Once  we  restrict  our  attention  to  x  in  a  compact 

£ 

K  c  ft,  however,  we  can  achieve  b  (0)  =  0  without  imposing  any 
further  assumptions.  The  following  argument  accomplishes  this: 
from  the  stability  of  0  with  respect  to  b°  as  in  d)  above  and  the  uniform 

convergence  of  bE  to  b^  one  can  deduce  that  (for  sufficiently 

£  £  £ 
small  £)  b  has  a  critical  point  C  such  that  5  0  as 

e  i  0.  Change  variables  to  y  =  x  -  £E.  The  new  coefficients 

ae(y)  =  ae(y+S£)  and  bE(y)  =  bE(y+£E)  satisfy  all  of  our 

"•*  £ 

assumptions  above  as  well  as  b  (0)  =  0.  The  only  difficulty  is 

£ 

that  the  domains  ft  -  S  are  e-dependent.  We  can  pass  to  a 

subdomain  ft'  so  that  y  €  ft'  implies  x  =  y  +  4>E  €  ft  and  a 

£ 

compact  subset  K'  c  q»  so  that  x  £  K  implies  y  =  x  -  C  £  K'  , 

for  sufficiently  small  e.  Applying  Theorem  1  to  ft',K'  we  get 

the  same  result  as  for  ft,K.  Taking  the  details  of  this  argument 

£ 

for  granted,  we  assume  in  the  following  that  b  (0)  =  0  and 
consequently  that  (2.2)  is  true. 

Here  is  our  main  theorem. 

Theorem  1:  Under  the  assumptions  described  above,  for  any  compact 
K  c  there  exists  6  >  0  and  £q  >  0  so  that  for  all  0  <  £  <  £ 


8 


sup  |uE(x)  -  u£(y)|  <  e  <S'£. 
x,y€K 

Roughly,  the  reasoning  behind  the  proof  is  that  for  small  e 

x  (t)  should,  with  high  probability,  follow  the  deterministic 

trajectory  x^(t)  into  the  vicinity  of  the  origin  before  making 

its  first  excursion  to  the  boundary  9ft.  A  precise  probabilistic 

estimate  along  these  lines  is  established  in  the  next  two  sections. 

To  apply  the  probabilistic  estimate,  we  need  to  know  a  modulus  of 

£ 

continuity  for  u  .  The  following  lemma  establishes  the  modulus 
that  we  need;  the  rescaling  argument  is  the  same  one  used  by 
Kamin  [ 5 ] . 

Lemma  1 :  Let  K  c  ft  be  compact.  Then  there  is  a  constant  C 
so  that  lVue(x)|  <  Ce’1^2  for  all  x  £  K  and  e  <  1. 

Proof : 

•  ^  “1/2  £ 

Make  the  change  of  variab^es^y  =  e  x.  Then  v  (y)  = 
ue(e^2y)  satisfies 

\  l  i (X)V  +  l  b£(y)V  =0  for  y  £  c'1/2ft. 
i,j  7rj  i  7  i 

The  coefficients  a£(y)  =  ae(e1/,2y),  be(y)  =  e"1//2be(e1//2y)  are 
Holder  continuous  with  respect  to  y  uniformly  in  e.  If  we  take  r  so  that 
Br(x)  =  { z :  |  x  -  z  |  <  r}  <=  ft  whenever  x  £  K,  then  Br(y)  <=  e~1//2ft 
whenever  y  e  e_1/2k  and  e  <  1.  We  can  apply  the  basic  Scahuder  interior 
estimate  ([3],  Theorem  6.2)  to  (y )  for  any  y  e  z'l^2Yi  to 
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conclude  that  |VvCy)|  <  C  for  all  y  €  K  for  some  constant  C. 
This  implys  the  lemma  after  changing  back  to  the  original 
variable  x. 
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III:  A  Prototype:  An  Ornstein-Uhlenbeck  Process 

Before  proving  the  estimate  on  hitting  probabilities 
of  the  next  section,  it  is  convenient  to  look  at  the  special  case 
of  an  Ornstein-Uhlenbeck  process  with  generator  as  in  (3.1)  below. 

The  proof  of  the  general  case  rests  on  a  comparison  with  the 
function  described  by  (3.4)  and  analyzed  below. 

In  1R^,  d  >  2,  suppose  that  a  >  0  is  a  constant  and  £E(t) 
is  a  diffusion  process  with  differential  generator 

(3.1)  yE[u](x)  =  j  Au(x)  -  ax-Vu(x). 

Let  T(r)  be  the  hitting  time  of  the  sphere  of  radius  r: 

(3.2)  x(r)  =  inf (t  >  0:  |CE(t)|  =  r). 

Take  a  fixed  R  >  0  and,  for  rQ  <  |x|  <  R,  define  the  hitting 

probability 

QvqM  =  Px[T(r0)  <  T(R);  T(r0)  <  °°1- 

What  we  show  is  that  there  exists  a  positive  constant  6^  >  0 
so  that 

e  -<5  /e  -6,/e  , 

(3.3)  Qr  (x)  >  1  -  e  1  whenever  r^  >  e  and  |x|  <y  R. 

In  words,  for  |se(0)|  <  ^  R  we  can  let  r^  I  0  exponentially 
with  e’1  and  at  the  same  time  have  the  probability  that 
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T(rQ)  <  T  (R)  converging  to  1  exponentially  fast.  To  prove  this 
we  calculate  Qe  (•)•  By  symmetry  Qe  depends  only  on  r  =  |x|. 

-  ro  e  ro 

Thus  Q  (x)  =  Q(r)  where  &  [Q(r)]  =  0  with  Q(r~)  =  1>  Q(R)  =  0. 
r0  U 

0  =  ^e[Q(|x|)]  =  f  Q"  (r)  +  [|  -  ar]  Q '  (r) 

or 

(3.4)  Q" (r)  +  [|  -  ^  r]Q'(r)  =  0;  Q(rQ)  =  1,  Q(R)  =  0. 


For  the  above  3  =  d  -  1,  but  we  will  carry  out  the  calculation 
for  arbitrary  positive  constants  a,B.  Solving  (3.4)  gives 


Q(r)  -  1  - 


a  2 

’r  -M  S  ds 
s  e 

r0 


a  2 

-Be  s  ds 
s  e 


If  rQ  <  r  <  R  ,  then 


„  a  2  ,  „ 
f  -0  e  s  (3  a  l  2 


R 


I 


ot  2 

R  -0^s 
s  ^ee  ds 

r0 


< 


[R  '  Irl 

J2„  i  R  e£  ^ 


=  e 


R) 


a  1  „2 

'  '  1  r0¥- 


R  3 


(R)' 


1  aR2 

If  log(rQ)  >  -  —  •  (-g^-) ,  then  the  preceding  is 


_  a  1  2 

V. 


<  e 
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Consequently,  if  is  slightly  less  than  the  minimum  of 

and  lIV  (slightly  less  so  as  to  absorb  the  R6),  then  for 

-^l/t  1 

rQ  >  e  and  r  <  j  R  we  have 

-6  /e 

(^•5)  Q(r)  >  1  -  e  for  sufficiently  small  e  >  o. 
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IV .  The  Hitting  Probabilities  in  the  General  Case 

Next,  we  prove  that  an  estimate  like  (3.  3)  holds  for  the 

hitting  probabilities  of  the  process  x6(t).  (t(r)  now  denotes 

the  time  of  first  contact  with  the  ball  of  radius  r  about  the 
£ 

origin  for  x  (•).) 


Theorem  2:  For  any  compact  K  c  fi  there  exists  6  >  0  so  that 

for  some  >  0  and  all  0  <  e  < 

[ t  ( r0 )  <  xfi]  -  1  '  e  <5/E  whenever  x  €  K  and  rQ  >  e  ^E. 

Proof : 

We  will  first  make  an  argument  for  an  appropriate  neighborhood 
of  the  origin.  The  key  to  the  proof  is  to  use  not  the  standard 
Euclidean  norm  |x|  but  a  different  symmetric  positive  definite 
quadratic  form.  By  hypothesis,  the  matrix  B  =  [3b? (0)/3Xj ]  is 
stable.  Lyapunov's  Theorem  on  matricies  implies  that  there  exists 
a  unique  symmetric  positive  definite  matrix  V  which  solves 

BTV  +  VB  =  -I. 

T  1  /  2  2 

Define  P(x)  =  [x  Vx]  '  .  For  f  €  C  QR) ,  a  computation  gives 
that 


(4.1)  ir£tf(p(x))]  -  ftp),  f  l  a?,p  p  *  £'(p)-[|  I  ,f.p  +  V  b.p  1 

i.j  7  *J  Vi  1 V 


VP(x)  ■  ~  ,  PXjX. 


V. . 
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The  idea  is  to  effect  a  comparison  of  each  of  the  terms  of 
(P)J  with  those  of  ^e[Q(r)]  in  (3.4).  First, 


v  „e  „  „  _  xTVaeVx 

i  aiipx.px.  — T - 

xi  xj  x  vx 


which  is  bounded  above  and  below  away  from  0  on  Q  -  (0).  More¬ 
over,  these  bounds  can  be  taken  to  be  independent  of  £  sufficiently 
small  since  ae  -*•  a^  uniformly.  Thus,  there  exists  a  constant 
A  >  0  so  that 


(4.2) 


A  *  <  I  a^.P  P  <  A  in  fi  -  {0},  all  sufficiently  small  £. 

ij  ^  Xj 


Secondly , 


(4.3) 


l  ae  p  =  I  7  ae  V  ♦  I  iLVajvx  <  C 
ij  xixj  p  ‘i  p  p2  -  p 


for  a  positive  constant  C  (again  uniformly  in  e  sufficiently 
small).  Thirdly, 


l  b£p  -  xT^  1*1  -  xVfxi  +  x!vLb£_-b0) 

4  1  X.  P(x)  P(X)  P 


Using  b°(x)  *  Bx  +  o(|x|)  and  b£  (x)  -  b°(x)  =  |x|o(l)  from 

(2. 2), we  find  that 


1  biPx.  "  ^[2-VBx  +  o(|X|)  +  0(1)] 

1  ft 


p'[-  \  +  0(|x|)  +  0(1)] 


■  g  miui  ;  . 


} 


15 

(The  o  ( |  x  | )  is  an  |x|  -*•  0  and  is  independent  of  e.  The  o(l) 

is  as  e  -*•  0  and  is  uniform  in  x.)  The  second  equality  is  a  con¬ 

sequence  of  our  choice  of  V.  It  follows  that  for  some  D,R  and 
all  positive, 

(4.4)  l  bfpx  <  -Dp (x)  if  P (x)  <  R  and  e  <  eQ. 

(Also  restrict  R  so  that  x  £  Q  whenever  P(x)  <  R.)  Take 
a  =  DA  ^ ,  $  =  AC  and  then  Q(')  as  in  (3.2).  Since  Q’  <  0, 
(4.1)-(4.4)  combined  imply  that,  for  e  <  eQ  and  P(x)  <  R. 

(4.5)  ^£[Q(P(x))]  >  j  A*1*  (QM(P)  +  Q’(P)[^  -  |ADP]}  =  0. 

Using  T (r)  for  the  hitting  time  of  the  set  {x:  P(x)  =  r)  by 

x  (•)»  (4.5)  implies,  either  via  the  maximum  principle  or  the 
fact  that  Q(P(xe(t)))  is  a  submartingale,  that 

P*[*(r0)  <  ?(R)]  >  Q(P(x)). 

If  Y  >  0  is  an  constant  so  that  Y  <  <  Y1,  then 

x  (Yrg)  >  x(rQ)  provided  |x£(0)|  >  Tq.  If  |x£(0)|  <  YR,  then  P(x£(0))  <  R  and 

x  (R)  <  xfi.  Consequently,  for  Tq  <  jx£(0)|  <  YR  we  have 

P^T(ro>  <  V  >  pxf?(Yro)  <  f(R)]  - 

This  is  trivially  true  also  if  |x£(0)|  <  r^.  The  calculation  of 

Section  2  now  implies  the  existence  of  >  0  so  that  for  all 
0  <  e  <  eg  and  |  x|  <  R 
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(4.6) 


-6^/e  -<5^/e 

Px[T(r0)  <  TnI  -  1  '  e  if  Vrn  >  e  1 


0 


The  last  step  of  the  proof  is  to  show  that  such  an  estimate 
remains  true  for  all  x  €  K.  By  the  strong  Markov  property 


pxerx(V  <  V  =  ex[p%(t(i  R))[I(ro)  <  T(-s R)  <  V 


*6l/e  e  y 

>  d-e  1  )-Px[t(}  R)  <  xfl]. 


It  is  sufficient  therefore  to  prove  that  for  some  6^  >  0  and 
all  x  €  K, 

e  Y 

(4.7)  Px[T(I  R)  <  V  >  1  '  e 

£ 

Let  <J>  (t;x),  c  >  0,  denote  the  solution  of  the  deterministic 
equation  <J>  *  ( t )  =  b  e  (4>  ( t ) )  with  $(0)  =  x.  Since  K  is 
contained  in  the  domain  of  attraction  of  the  stable  point  0, 
there  exist  T,n  >0  so  that  if  x  €  K  then 

y  €  U  whenever  |y  -  <P°(s;x)|  <  2n  for  some  0  <  s  <  T, 

and 

| y |  <  |  R  whenever  |y  -  $°(T;x)|  <  2n. 

As  e  -*•  0,  be  converges  to  b^  uniformly  in  q  and  consequently 

$e(t;x)  converges  to  $°(t;x)  uniformly  for  t  £  [ 0 ,  T]  and  x  €  K. 
Therefore,  if  e  is  sufficiently  small  and  x  £  K  it  will  be  true 
that 
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y  €  fi  whenever  |y  -  4> C ( s ; pc )  |  <  n  for  some  0  <  s  <  T, 

and 

|y|  <  j  R  whenever  | y  -  $  (T;x)|  <  n. 

For  such  £  and  x  6  K, 

P£[x0  <  T (t  r)3  <  p!t  SUP  l<fe(s,x)  -  Xe(s)|  >  n]. 

X  w  -  J  X  0<s<T 

Define 

0£(t)  =  x£(t)  -  x  -  (  be(x£(s))ds 

J  0 

/  t  £  £  a 

/c  j  o  (x  (s ) ) ds  if  x  is  obtained  from  an  Ito  equation 
Gronwall's  inequality  implies  that 


sup  |4>E(s;x)  -  xE(s)|  <  eMT  •  sup  |6E(s)| 

[0 ,T]  ‘  0<s<T 

where  M  is  the  Lipschitz  constant  for  the  b E  C  * )  (uniform  in  £) . 
It  is  a  standard  argument,  using  exponential  martingales,  that 

P£[  sup  |9£(s)|  >  Xj  <  (2d)exp[-2-;Lf] 

x  j0>T]  Zd£Ai 

T  2 

where  x  acx  <  Aj|x|  j  ;  see  [11],  equation  (2.1),  pg.  87  and  proof  for 
instance.  Combining  these  facts,  for  all  x  €  K  and  £ 
sufficiently  small, 


L 
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E  y  1  n2e ' 2MT 

pxITf!  S  R)1  ?  C2d)exp[-  i 

This  shows  (4.7)  and  completes  the  proof. 

Theorem  1  is  now  simple. 

Proof  (of  Theorem  1): 

Take  any  x  6  K  and  set  r^  =  e  '  (6  >  0  as  in  Theorem  2). 

uE(x)  -  E^[ue(xe(xfi(r0)));  T(rQ)  <  xfl]  +  E^[ f £ (x£(Tfi) ) ;  <  t(r0 

Therefore , 

uE(x)  -  ue(0)  =  E£[ue(x£(Tfi(r0)))  -  uE(0);  x(r0)  <  xfl] 

+  Exe[f£(xe(Tfl))-  ue(0);  <  x  (r  Q )  ] 

(ue(x)  -  U£(0)|  <  sup  |ue(y)  -  U£(0)|  + 

|y|<r0 

2  s3uPt£el,px[Tfi  *  nvi 

<  sup|Vue|*rn  +  2  sup)  f £)  e"  5//£’ 

K  u  an 

<  ce^V^+Z  suP|fe)e-5/£. 

an 


The  theorem  now  follows  (with  a  new  slightly  smaller  <5) . 
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V :  Concluding  Remarks 

We  have  two  simple  observations  to  make  in  closing.  The  first 

is  regarding  the  case  in  which  ft  contains  several  critical  points 

of  (1.1).  This  has  been  discussed  in  the  literature:  [9],  [12]. 

£  £ 

Both  u  -*■  constant  and  u  -»•  a  piecewise  constant  function  are 
possibilities  now,  depending  on  the  Ventcel-Freidlin  variational 
distances  between  the  critical  points  and  3ft.  If  x*  is  an 
asymptotically  stable  critical  point  (replacing  the  origin  in  d) 
of  section  2)  and  ft*  c  ft  is  its  domain  of  attraction,  then  by 
taking  fe  =  u£  on  3ft*  we  can  apply  Theorem  1  to  see  that 
leveling  takes  place  exponentially  fast  in  each  such  domain  of 
attraction. 

Finally,  we  observe  that  the  specification  of  boundary  data 

£ 

f  is  actually  superfluous.  All  that  matters  in  the  proof  is  the 
availabili  y  of  a  bound  in  e  for  the  uE.  Theorem  1  could  be 
formulated  as  follows : 

for  K  c  ft  compact  there  exist  6  >  0  and  eQ  >  0  so  that 
£ 

whenever  ¥  [u]  =  0  in  ft  and  e  <  eQ , 

(6.1)  sup  |u(x)  -  u(y)  |  <  e‘6/,£sup|u|  . 

x,y€K  "  ft 

Define  the  exit  measures  on  the  Borel  subsets  of  3ft  by 


*>)  «  P>E(V  €  B]. 


The  strong  maximum  principle  implies  that 


and 


are  mutually 


in*  it— 
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